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We propose PALPS, a Process Algebra with Locations for Population Systems. PALPS allows us to 
produce spatially-explicit, individual-based models and to reason about their behavior. Our calculus 
has two levels: at the first level we may define the behavior of an individual of a population while, 
at the second level, we may specify a system as the collection of individuals of various species 
located in space, moving through their life cycle while changing their location, if they so wish, and 
interacting with each other in various ways such as preying on each other. Furthermore, we propose 
a probabilistic temporal logic for reasoning about the behavior of PALPS processes. We illustrate 
our framework via models of dispersal in metapopulations. 

1 Introduction 

During the last decade we have witnessed an increasing trend towards the use of formal frameworks for 
reasoning about biological as well as ecological systems including process algebras |[3Tl[30l l9l l24l[T4l . 
Membrane Systems ||29l[TTI and cellular automata ifTTl . Process algebras, first proposed in ||25ll20l to aid 
the understanding and reasoning about communication and concurrency, provide a number of features 
that make them suitable for capturing biological processes. In particular, process algebras are especially 
suited towards the so-called "individual-based" approach of modeling populations, as they enable one to 
describe the evolution of each individual of the population as a process and, subsequently, to compose a 
set of individuals (as well as their environment) into a complete ecological system. Features such as time, 
probability and stochastic behavior, which have been extensively studied within the context of process 
algebras, can be exploited to provide more accurate models, while associated analysis tools can be used 
to analyze and predict their behavior. 

In this work, our aim is to introduce a process-algebraic framework to enable spatially-explicit mod- 
eling of ecological systems. Such modeling |[T5l l31 has been of special interest to conservation scientists 
and practitioners who have employed it in order to predict how species will respond to specific manage- 
ment schemes and guide the selection of reservation sites and reintroduction efforts, e.g. |[T9ll28l . The 
use of spatially-explicit, individual-based modeling requires the description of the environment and the 
individuals residing in it, including a description of each individual's interaction with other individuals 
as well as with the environment. As far as the environment is concerned, these models typically involve 
the use of patches or a lattice to represent the habitat. Individuals are then placed on specific locations 
of the modeled landscape and their behavior, including events such as birth, mortality, and dispersal, is 
simulated at the individual or the population level and analyzed. 

In order to capture this type of behavior our process algebra, PALPS, associates processes with 
information about their location and their species. The habitat is defined as a graph consisting of a set 
of locations and a neighborhood relation. Movement of located processes is then modeled as the change 
in the location of a process, with the restriction that the originating and the destination locations are 
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neighboring locations. In addition to moving between locations, located processes may communicate 
with each other by exchanging messages upon channels. Communication may take place only between 
processes which reside at the same location while special channels allow processes to engage in preying 
and reproduction. Furthermore, PALPS may model probabilistic events, with the aid of a probabilistic 
choice operator, and uses a discrete treatment of time. Finally, in PALPS, each location may be associated 
with a set of attributes capturing relevant information such as the capacity or the quality of the location. 
These attributes form the basis of a set of expressions that refer to the state of the environment and are 
employed within models to enable the enunciation of location-dependent behavior. 

The operational semantics of our calculus is given in terms of a labeled transition system on which 
we may check properties expressed in an instantiation of the PCTL temporal logic. We illustrate the 
expressiveness of PALPS by constructing spatially-explicit individual-based models for metapopulation 
dispersal. 

There exists a variety of previous proposals which introduce locations or compartments into formal 
frameworks, e.g. (21 HQl [T3l [27l [22l |4l [7J , while work has been carried out to employ these frameworks 
for modeling and analyzing population systems [5]. PALPS departs from these works in that it is the 
first process-algebraic framework developed specifically for reasoning about ecological models as well 
as in its treatment of a state and its capability of expressing state-dependent behavior. In particular, 
it can be considered as an extension of WSCCS of 0T1 with locations and location attributes, while it 
shares a similar treatment of locations with process algebras developed for reasoning about mobile ad hoc 
networks, e.g. ll23l[T8l . As such, PALPS considers a two-dimensional space where locations and their 
interconnections are modeled as a graph upon which individuals may move as computation proceeds. 
The main feature that distinguishes PALPS from existing formal frameworks is the fact that it associates 
locations with a set of attributes that model special characteristics of locations which may be of interest 
when modeling a system and the ability to express behavior of individuals that is conditional on the 
values of these attributes. Examples of attributes that can be observed by individuals is the number of 
individuals a location can support as well the current number of individuals present at a location. 

In the remainder of the paper we present the syntax and the semantics of PALPS in Section 12 while 
in Section[3]we provide models of metapopulation dispersal. In Section|4]we conclude with a discussion 
on future work. 

2 The Process Calculus 

In our calculus, PALPS (Process Algebra with Locations for Population Systems), we consider a system 
as a set of individuals operating in space, each possessing a species and a location identifier. Movement 
in the calculus is modeled via a specialized action whose effect is to change the location of an individual, 
with the restriction that the originating and the destination locations are neighboring locations. The 
notion of neighborhood is implemented via a relation Nb where (£,£') G Nb exactly when locations £ 
and £' are neighbors. We also use Nb as a function and write Nb(^) for the set of all neighbors of £. 

2.1 The Syntax 

We continue to formalize the syntax of PALPS. We begin by describing the basic entities of the calculus. 

• We assume a set of special labels S corresponding to the species under consideration, ranged over 
by s, s'. 
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• Furthermore, we assume a set of channels Ch, ranged over by lower-case strings. This set contains 
the special channels rep s and prey s , s 6 S, which are channels used to model reproduction of 
species s and preying on species s. 

• Finally, we assume a set of locations Loc ranged over by £, £'. Locations can be associated with 
a set of attributes that model special characteristics of locations of interest within a system. We 
write Y f° r attributes and f° r the value of attribute y at location £. 

Our calculus also employs two sets of expressions: logical expressions ranged over by e and arith- 
metic expressions, ranged over by w. One of our main aims being to facilitate reasoning about spatially- 
dependent behavior, these expressions are intended to capture environmental (location-relevant) situa- 
tions which may affect the behavior of individuals. Expressions e and w, are constructed as follows: 

e ::= true | -ie | e\ Ae 2 \ w co c 

w ::= c | y@t \ s@t \ @t | op^w) | op 2 (wi,w 2 ) 

where c is a real number, cxig {=,<,>} and £* G LocU {myloc}. Let us informally consider the intro- 
duced expressions. To begin with, logical expressions are built using the propositional calculus connec- 
tives as well as comparisons between an arithmetic expression w and a constant c, i,e. w \x\ c. Moving 
on to arithmetic expressions, these include three special expressions interpreted as follows: Expression 
y@£* is equal to the value of attribute y at location £*. Expression (s@£*) is equal to the number of 
individuals of species s at location £* and expression @£* denotes the total number of individuals of all 
species at location £*. As specified above, 1* can be an arbitrary location or the special location myloc. 
This label is employed to bestow individuals the ability to express conditions on the status of their current 
location no matter where that might be as computation proceeds. Specifically, myloc refers to the actual 
location of the individual in which the expression appears and it is instantiated to this location when the 
condition needs to be evaluated (see rule (Cond) in Table [3]). 

Thus, arithmetic expressions are the set of all expressions formed by arbitrary constants c, quantities 
Y@£*, s@£*, @£*a.nd the usual unary and binary arithmetic operations (opj and op 2 ) on the real num- 
bers. Logical expressions and arithmetic expressions are evaluated within a system environment. The 
precise definition of the evaluation function is postponed to Tables Q] and |2] 

We may now move on to the syntax of PALPS which is given at three levels: (1) the individual level, 
ranged over by P, (2) the species level, ranged over by R, and (3) the system level, ranged over by S. 
Their syntax is defined via the following BNF's: 

P ::= | T].P | J^Wi-.P | cond (e 1 >Pi,...,e n >P n ) | C 
iel 

R ::= \rep.P 

S ::= | | 7?:[[s]] | Si\S 2 \ S\L 

where a € Ch, L C Ch, C ranges over a set of process constants ^, each with an associated definition of 

def 

the form C = P, where the node P may contain occurrences of C, as well as other constants, and 

r] ::= a \ a \ go £ \ . 

Beginning with the individual level P, process represents the inactive individual, that is, an indi- 
vidual who has ceased to exist, tj .P describes the individual who first engages in activity r\ and then 
behaves as P. Activity r\ can be an input action on a channel a, written simply as a, a complementary 
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output action on a channel a, written as a, a movement action with destination £, got, or the time-passing 
action, y/ . Actions of the form a, and a, a G Ch, are used to model arbitrary activities performed by an 
individual e.g. eating, preying, observing the environment as well as reproduction. Thus, for example, 
the actions prey s and prey s are executed, respectively, by a prey of population s and a predator who is 
preying on individuals of population s. The tick action y/ measures a tick on a global clock and is used 
to separate the phases/rounds of an individual's behavior. Essentially, the intention is that in any given 
time unit all individuals perform their available actions possibly synchronizing as necessary until they 
synchronize on their next y/ action and proceed to their next round. 

Y<iei w i '■ Pi represents the probabilistic choice between processes P„ i 6 /. Each alternative is asso- 
ciated with a probability of appearance, which is the value to which the expression w; evaluates. The 
conditional process cond (e\ >P\, . . . , e n >P n ) presents the conditional choice between a set of processes: 
it behaves as Pi, where i is the smallest integer for which e, evaluates to true. Finally, process constants 
provide a mechanism for including recursion in the calculus. 

Moving on to the action of reproduction, to capture the creation of new individuals, we employ the 
special species processes R. R, defined as Irep.P, are replicated processes which may continuously re- 
ceive input through channel rep and create new instances of process P, where P is a new individual of 
species R. Such inputs will be provided by individuals in the phase of reproduction via the complemen- 
tary action rep. 

Finally, population systems are built by composing in parallel located individuals, P:[[s,£]], where s 
and t are the species and the location of the individual, and species P:[[s]], where s is the name of the 
species. Finally, S\L models the restriction of the use of channels in set L within 5. 

As an example, we consider the model described in JH where a set of individuals live on an n x n 
lattice of resource sites and go through phases of reproduction and dispersal. Specifically, the studied 
model considers a population where individuals disperse in space while competing for a location site 
during their reproduction phase. They produce an offspring only if they have exclusive use of a location. 
After reproduction the offspring disperse and continue indefinitely with the same behavior. In PALPS, 

def 

we may model the described species s as R = Irep.P, where 

def 1 

P = - : go £.y/. cond (s@myloc = 1 >P\\ true> y/.P) 

f.£Neigh(my\oc) 

def 

Pi = p : rep.y/.P\ + {l—p): rep.rep. y/.P\ 

We point out that the conditional construct allows us to determine the exclusive use of a location by an 
individual. The special label myloc is used to illustrate that the location of interest is the actual location 
of an individual once the individual is placed in a context within a system definition. Furthermore, note 
that Pi models the probabilistic production of one or two offsprings of the species. During the dispersal 
phase, an individual moves to a neighboring location which is chosen probabilistically among the four 
neighboring locations on the lattice of the individual. Then a system containing of two individuals at a 
location I and one in location £' can be modeled as 

^em d ^(P:[[As]]|P:[[£, S ]]|P:r, S ]]|(!r^.P):[[ S ]])\{r^}. 

def 

To model a competing species s' which preys on s, we may define the process R f = \rep f .Q, where 
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cond (s@myloc > 1 r> prey s .y/.Qi, true > y/.Qz) 
rep'.^.Q 

cond (s@myloc > 1 r> prey s . y/ '.<2i, true [> -y/.O) 

An individual of this species looks for a prey. If it succeeds in locating one, then it produces an offspring. 
If it fails for two consecutive time units it dies. 

2.2 The Semantics 

The semantics of PALPS is defined in terms of a structural operational semantics given at the level 
of configurations of the form (E,S), where E is an environment and S is a population system. The 
environment E is an entity which captures how the various locations of the system are populated. More 
precisely, E C Loc x S x N, where each pair £ and s is represented in E at most once and where (£,s,m) 6 
E denotes the existence of m individuals of species s at location £. The environment E plays a central role 
in defining the semantics of the calculus and, in particular, for evaluating expressions. The satisfaction 
relation for logical expressions |= is defined inductively on the structure of a logical expression as shown 
in Tabled] 

Table 1 : The satisfaction relation for logical expressions 



Zs^true 


always 




E |= 


if and only if 




E |= e\ Ae 2 


if and only if 


E \= e\ t\E = e2 


E\=wx\e 


if and only if 


val (E,w) tx e 



Q 

Qi 



def 



def 



def 



The relation |= is straightforward and depends on the evaluation function for arithmetic expressions 
val (E,w) defined in Table [2] 



Table 2: The evaluation relation for arithmetic expressions 



val (E,c) 
val {E,yf@e) 
val (E,s@£) 
val (E,@e) 
val (£,op!(w)) 
val (£ , ,op 2 (wi, wi)) 



c 

We 

num(£,f,s) 

num'(£,l) 

opi(val (E,w)) 

op 2 (val (E,wi),va\ (E,W2)) 



The auxiliary functions num(£,l,s) and num'(E,£) compute the number of individuals at location £ 
in environment E of a specific species s (num(£,£,s)) or for all species (num'(£,f)) and are defined by 
num(£,l,s) = n where (£,s,n) G E and num'(E,£) = Y,ses num(E,s,£). 

Before we proceed to the semantics we define some additional operations on environments that we 
will use in the sequel: 



M. Antonaki and A. Philippou 



19 



Definition 1. Consider environment E location £ and species s. 

• E® (s,£) increases the count of individuals of species s at location £ in environment E by 1; 



• EQ (s,£) decreases the count of individuals of species s at location I in environment E by 1: 



We may now define the semantics of PALPS, presented in Tables [3] and [4] and given in terms of two 
transition relations, the nondeterministic relation — > n and the probabilistic relation — > p . A transition of 

the form (E,S) (E',S') signifies that configuration (E,S) may execute action p and become (E',S f ) 
whereas a transition of the form (E,S) —*p (E',S') signifies that configuration (E,S) may evolve into 
configuration (E',S') with probability w. Whenever the type of the transition is irrelevant to the context, 

we write (E,S) -A (E',S') to denote that either (E,S) (E',S') or (E,S) (E',S'). Action p 
appearing in the nondeterministic relation may have one of the following forms: 

• a@£ and a@£ denote the execution of actions a and a respectively at location £. 

• z denotes the internal action. This may arise when two complementary actions take place at the 
same location or when a move or a prey action take place. We are not interested in the precise 
location of internal actions, thus, this information is not included. 

• sj denotes the time passing action. 

The rules of Table [3]prescribe the semantics of located individuals in isolation. The first four axioms 
define nondeterministic transitions, the fifth axiom defines a probabilistic transition, and the last two 
rules refer to both the nondeterministic and the probabilistic case. All rules are concerned with the 
evolution of the individual in question and the effect of this evolution to the system's environment. A 
key issue in the enunciation of the rules is to preserve the compatibility of P and E as transitions are 
executed. We consider each of the rules separately. Axiom (Tick) specifies that a ^-prefixed process 
will execute the time consuming action yj and then proceed as P. The state of the new environment 
depends on the state of P: if P = then the individual has terminated its computation and, therefore, it 
is removed from E (see the definition of E p ) whereas, if P / then, obviously, E remains unchanged. 
Axiom (Act) specifies that tj .P executes action r\ @£ and evolves to P. Note that the action is decorated 
by the location of the individual executing the transition to enable synchronization of the action with 
complementary actions taking place at the same location (see rule (Par2), TableHJ). This axiom excludes 
the case of r\ = go £ which is treated separately in the next axiom. Specifically, according to Axiom 
(Go), an individual may change its location. This gives rise to action % and has the expected effect 
on the environment E. Moving on to Axiom (Prey), this describes that any individual can become 
the victim of a preying action. This may happen at any point during the lifetime of the individual 
giving rise to the action prey s @£ and causing the individual to terminate with the appropriate changes 
to the state of the environment. Rule (PSum) expresses the semantics of probabilistic choice: once the 
probability expressions are evaluated within the environment, the probabilistic action is taken leading to 
the appropriate continuation: if the resulting state of the individual, namely Pi, is equal to 0, then the 




E' U { (£, s, m + 1 ) } ifE = E' U { (£, s, m)} for some m 
EL){(£,s, 1)} otherwise 
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Table 3: Transition rules for individuals 



(Tick) (E,^.P:l Sl £]]) (£V:[M]) 

(Act) {E,r,.P:ls,q^i n (E p ,P:ls,q) r]^go£' 

(Go) {E,got'.P:ls,q) ((£eM))©M'),P:[M']]) G Nb 

(Prey) (£,P:M pr ^ n (Ee(s,£),0:^A) 

(PSum) (£,L e /w,- : ^[MD val ^p% (£V;:[[M]]) 



(Const) (^MD^'^MII) c^P:M 

(Cond) (gv^jMD {E' ^-.[[s/WlE^e^E HjjlJ < i 
(£,cond (ei>Pi,...,e n >P n ))^{E',P> i :ls,£'l) 



£ otherwise 



individual is removed from the environment E. Note that we write w\.l for the expression w with all 
occurrences of myloc substituted by location I: w\.£ = w[£/myloc]. Next (Const) express the semantics 
of process constants in the expected way. Finally, rule (Cond) stipulates that a conditional process may 
perform an action of continuation Pi assuming that ej\,£ evaluates to true and all ej\ r £, j < i evaluate to 
false. Similarly to w\,£, e\.£ is the expression e with all occurrences of myloc substituted by location I. 

We may now move on to Table|4]which defines the semantics of system-level operators. The first rule 
defines the semantics for the replication operator, the next five rules define the semantics of the parallel 
composition operator, and the last rule deals with the restriction relation. 

Thus, according to axiom (Rep), a species process may execute action rep s @£ for any location £ and 
create a new individual P of species s at location £. Next, rules (Pari) - (Par4) specify how the actions 
of the components of a parallel composition may be combined. Note that the symmetric versions of 
these rules are omitted. According to (Pari), if a component may execute a nondeterministic transition 
and no probabilistic transition is enabled by the other component (denoted by (E,S2)-/-^- p ), then the 
transition may take place. If the parallel components may execute complementary actions, then they 
may synchronize with each other producing action T (rule (Par2)). If both components may execute 
probabilistic transitions then they may proceed together with probability the product of the two distinct 
probabilities (rule (Par3)) and, finally, if exactly one of them enables a probabilistic transition then this 
transition takes precedence over any nondeterministic transitions of the other component (rule (Par4)). 
Note that in case that the components proceed simultaneously then the environment of the resulting 
configuration should take into account the changes applied in both of the constituent transitions (rules 
(Par2) and (Par4). This is implemented by E® (E\,E2) as follows: 



£<g) (Ei,E 2 ) = {(£,s,m + ii +/ 2 ) | (£,s,m) G E, {£,s,m + i\) G E\, (£,s,m + i2) £ E2,h,i2 6 
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Table 4: Transition rules for systems 



R=\re Ps .P:{{s]\,£eLoc 
{E,R) re ^' n (E . is./)./': s.l R) 

(£,5i)-^„ (E',S[),(E,S 2 )^ P 
(£,Si|S 2 ) A„(E',Si|S 2 ) 

(£,SQ^i„ (£ 1 ,y i ) ; (£,5 2 ) a g|„(£ 2 ,5 2 ) 
(£,Si|S 2 ) -^ n (£®(£i,£ 2 ),5i|5 2 ) 

(£,SQ A, (g^Cg^) -^> P (£2,^2) 
(E.SilSa) 1 ^, (£®(£i,£ 2 ),5;|5 2 ) 

(£,5i|5 2 )^(£',5i|5 2 ) 
(£,50 A, (£i,S / 1 ),(£,S 2 ) A, (£ 2; 5 2 ) 

(£,5) -A (E',S'),a<£ {a@e,a@e\g € L} 
(£,S\L) -A (E',!f)\L 



Next, rule (Time) defines that parallel processes must synchronize on sj actions, thus allowing one tick 
of time to pass and all processes to proceed to their next round. Finally, rule (Res) defines the semantics 
of the restriction operator in the usual way. 

Based on this machinery, the semantics of a system S is obtained by applying the semantical rules 
to the initial configuration. The initial configuration, (E,S), is such that (£,s,m) 6 E if and only if S 
contains exactly m individuals of species s located at I. In general, we say that E is compatible with 
S whenever (£,s,m) £ E if any only if S contains exactly m individuals of species s located at I. It is 
possible to prove the following lemma by structural induction on S fill. 

Lemma 1. Whenever (E,S) — > (E',S') and E is compatible with S, then E' is also compatible with S'. 
2.3 Model Checking PALPS 

Model-checking of PALPS processes may be implemented via an instantiation of the PCTL logic (H. 
The instantiation involves the adoption of PALPS logical expressions as the atomic propositions of the 
logic. Specifically, the syntax of the PCTL instantiation that we consider, is given by the following 
grammar where <1> and <p range over PCTL state and path formulas, respectively, p £ [0, 1] and k £ N. 

<E> := true \ e | -i3> | <J>A<J>' | P^pl^] 
:= X<£ I <£U*<!> I <J>iUO 

In the syntax above, we distinguish between state formulas <1> and path formulas 0, which are eval- 
uated over states and paths, respectively. A state formula is built over PALPS logical expressions and 
the construct P Mp [0]. Intuitively, a configuration s satisfies property P M p[0] if for any possible execution 



(Rep) 

(Pari) 

(Par2) 

(Par3) 

(Par4) 

(Time) 

(Res) 
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beginning at the configuration, the probability of taking a path that satisfies the path formula satisfies 
the condition oo p. Path formulas include the X (next), (bounded until) and U (until) operators, which 
are standard in temporal logics. Intuitively, X<J> is satisfied in a path if the next state satisfies path formula 
<J>, <J>i U k <t>2 is satisfied in a path if <J>i is satisfied continuously on the path until <I>2 becomes true within 
k time units (where time units are measured by y/ events in PALPS) and is satisfied if <I>2 is 

satisfied at some point in the future and <I>i holds up until then. 

For example, consider a population s in danger of extinction. A property that one might want to check 
for such a population is that the probability of extinction of the population in the next ten years is less than 
a certain threshold p e . This can be expressed in PCTL by the property P < Pe [true\} w Y.teLoc s @^- = 0]- 
Alternatively, one might express that a certain central location t will be reinhabited with at least some 
probability p r by: = 0—7- P> Pr [true\J(s@£ > 0)]. Similarly, it would be possible to study the 
relation within a model between the size of the initial population and the probability of extinction of 
the population, by checking properties of the form > m — > P> Pe [true\}{s@l = 0)] or explore the 
dynamics between two (or more) competing populations s and s' by, for example, expressing that within 
the next 20 years with some high probability, members of the population s will outnumber the members 
of population s': P>p[trueUC£ £eIj0c s' @£ < E^ eL oc s@ ^)]- 

The semantics of PCTL are defined over Markov Decision Processes (MDPs), a type of transition 
systems that combine probabilistic and nondeterministic behavior. It is not difficult to see that the opera- 
tional semantics of PALPS gives rise to transition systems that can easily be translated to MDPs HI. For 
the details of the semantics and the model checking algorithm we refer the reader to |[T6ll . 

As a final note we observe that in order to check the satisfaction of PCTL properties by PALPS pro- 
cesses it is sufficient to restrict our attention to the E component of each configuration (E,S). This is due 
to the fact that E is the only information required in order decide the satisfaction of logical expressions 
by configurations (see Tables Q] and [2]). 

3 Examples 

During the last few decades, the theory of metapopulations has been an active field of research in Ecol- 
ogy and it has been extensively studied by conservation scientists and landscape ecologists to analyze 
the behavior of interacting populations and to determine how the topology of fragmented habitats may 
influence various aspects of these systems such as local and global population persistence and species 
evolution. The notion of a metapopulation refers to a group of distinct populations of the same species 
residing on a fragmented habitat or, a so-called set of patches, and cycle in relative independence through 
their life cycle while interacting with other populations and colonizing previously unoccupied locations 
through dispersal. It has been observed that while populations of a metapopulation may go extinct 
as a consequence of demographic stochasticity, the metapopulation as a whole is often stable because 
immigrants from another population are likely to re-colonize habitat which has been left open by the 
extinction of another population or because immigration to a small population may rescue that popula- 
tion from extinction. Indeed the process of dispersal is of vital importance in metapopulations. It affects 
the long-term persistence of populations, the coexistence of species and genetic differentiation between 
subpopulations and understanding this process is essential for obtaining a good understanding of the be- 
havior of metapopulations. The evolution of dispersal has received much attention by scientists and it 
has been studied in connection to various parameters such as the connectivity of the habitat on which a 
metapopulation exists, patch quality and local dynamics. 

In this section, we describe two examples relating to metapopulation dispersal through which we 
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Figure 1 : The sequence of events in the lifetime of a dispersing species 



illustrate how our calculus can be used to construct models of this phenomenon. 

Example 1. The first example we consider is motivated by the spatially-explicit, individual-based model 
of IT321 . In this work the authors construct a fairly simple model of metapopulation dispersal which 
departs from previous works in that, unlike previous models of metapopulation dispersal which tended 
to be deterministic and at the level of population densities, the model constructed is both stochastic and 
individual-based. 

According to this study, a set of genotypes co-exist within a habitat which differ only in their propen- 
sity to disperse. The metapopulation is composed of n x n subpopulations inhabiting a set of patches 
arranged on a square lattice with cyclic boundaries, so that individuals leaving the "top" or "right-side" 
of the world reappear on the "bottom" or "left-side" respectively and vice versa. Each patch is associated 
with a so-called patch quality related to the capacity of the patch. The behavior of an individual of the 
genotypes under study is illustrated diagrammatically in Figure Q] According to this model, an adult 
individual initially produces X offsprings. Subsequently, a phase of competition takes place between the 
juveniles of the population of which a fraction survives. Each surviving offspring may disperse accord- 
ing to a probability of dispersal distinct to its genotype. In case it disperses, the neighboring patch it 
moves to is selected with equal probability among all neighbors. This sequence of events in the behavior 
of an individual is presented diagrammatically in Figure [TJ We point out that the percentage of offspring 
surviving juvenile competition at patch I is given by ye = (1 + CLt • Ni)^, where CLg is the measure of the 
patch quality, N& is the number of individuals residing at patch I and j3 is a constant that relates to the 
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degree of competition. 

This metapopulation can be modeled in PALPS as follows. We consider the set of of locations (/, j), 
1 < i,j < n, where two locations and (k,l) are neighbors if they are adjacent on the grid. Finally, let 
us consider the location attribute Cty as a measure of the quality of the patch at t. Then, genotype i with 
some constant probability of dispersal pi and A = 3 can be defined as the species process Rj =\rep\.Ju 
where 

repj.fepj.repi.O Adult Individual 

qi : Si + (1 — qi) : Juvenile 
Pi : Dj + (1 — pi) : y/Aj Surviving Juvenile 

LeeNei g h( my \oc) \-got V A i Dispersing Juvenile 

and qi the probability of survival of juvenile competition is given by q i =(l + a r @£)P. Then a system 
can be modeled as the composition of the various genotypes as well as the individuals of the initial 
population under study: 

System = [(RrM I ---RkM I EI ^Pi' 1 ! I ■ ■ -)\{repu ■ ■ -re Pk }. 

\<i<m\ 

Analysis in this model may focus on the effect that the dispersal rates, the degree of competition 
and/or patch quality may have on the degree of population dispersals. 

Example 2. As another more complex example, let us consider a model of wood thrush dispersal, 
initially proposed in 11331 and expanded upon in [26]. This model considers three types of birds: adult 
breeders, adult floaters, and juveniles which are birds in their first year of life. According to this model, 
adult breeders produce an offspring at a rate dictated by various system parameters such as clutch size, 
nest predation and paratisism rates which we denote as ry. Following reproduction, each individual 
has a probability of dying before the next time step which is higher in juveniles and adult floaters in 
comparison to adult breeders. We write qj,, qj and qf for the mortality rates of breeders, juveniles and 
floaters, respectively. If following mortality a habitat patch has more birds than its capacity allows, 
then dispersal will occur according to a probability determined by the size of the patch and the distance 
between neighboring patches. This probability is higher in floaters and juveniles in comparison to adult 
breeders who exhibit a high site fidelity. We write pb, Pj and pj for the dispersion rates of breeders, 
juveniles and floaters, respectively. If a bird reaches a patch with available capacity then it will settle. 
If not, then it will either attempt to disperse to another patch or it will become a floater depending 
on whether it has reached its maximum number of dispersal events. Once dispersal has occurred, the 
juveniles become adults and the model begins another cycle. This sequence of events in the behavior of 
the populations is presented diagrammatically in Figure [2] 

This metapopulation can be modeled in PALPS as follows. We consider the set of of locations and an 
associated predefined neighbor function as well as a distance function that may be instantiated according 
to the modeler's preference to capture Euclidean distance or some other function of interest [26] . We 
also assume the existence of a set of probabilities {pij}ijehoc where pij represents the probability of 
dispersal from patch i to patch j. Finally, we introduce the location attribute ce as a measure of the 
capacity of patch I. Then, wood thrush species can be modeled by the process R = \rep.Juv, where the 
behavior of a juvenile individual 7, is described by the following equations: 



A, 
Ji 
Si 
D, 



def 
def 
def 
def 
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Breeders 



Floaters 



Reproduction 



Juveniles 



Breeders 



No need for dispersal 
Surviving Juveniles 
And Breeders 



Non-dispersing 
juveniles and breeders 



No need for dispersal - 
Surviving Juveniles 
And Breeders 



Non-dispersing 
juveniles and breeders 



No need for dispersal ■ 
Surviving Juveniles 
And Breeders 



Mortality 



Floaters 



Surviving Breeders 
Juveniles and Floaters 



^No need for dispersal - 
Surviving Floaters 



Need for dispersal 



Dispersing 
juveniles, breeders 
and floaters 



^Non-dispersing 
floaters 



Dispersal 



No need for dispersal - 
Surviving Floaters 



Need for dispersal 



Dispersing 
juveniles, breeders 
and floaters 



"Non-dispersing 
floaters 



Dispersal 



Need for dispersal 



No need for dispersal - 
Surviving Floaters 



Breeders' 



Juveniles and Floaters 



Figure 2: A cycle in the lifetime of the metapopulation 
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Juv 


def 


qj : JC +(l-qj) : 


Juvenile survival 


JC 


def 


cond (@myloc > c my | oc [>/D ,true[> y/ AB) 


Check patch capacity 


JD 


def 


PjlJAi + il-p^'.y/AB 


Decide whether to disperse 


J A x 


def 


Y,i€Neigh(my\oc) Pmyloci : goi.JC\ 


Dispersal attempt 1 


JC\ 


oeT 


cond (@myloc > c my | oc > JD\ , true > y/ AB) 


Check patch capacity 


J L)\ 


def 


Pj . JA2 + (I — Pj) . V .AB 


Decide whether to disperse 


T A 

JA2 


def 


S~ n . ,,,, (I if 

LteNeigh(my\oc) Pmy\oc,i ■ gOl.JL,2 


Dispersal attempt 2 




def 


COnu y V"/ my IOC ,> ^-myloc ^ V ' trlic L> -y ./i/J J 


oecome noaier or aoun 


AB 


def 


r b :7bi.BS+(l-rb):BS 


Breeder reproduction 


BS 


def 


q b :BC + (l-q b ):0 


Breeder survival 


BC 


def 


cond (@myloc > c my | oc >BD ,true > y/AB) 


Check patch capacity 


BD 


def 


Pb-.BAx + il-p^-.y/AB 


Decide whether to disperse 


BA\ 


def 


Y,leNeigh{my\oc) Pmyloci : goi.BC\ 


Dispersal attempt 1 


BC X 


def 


cond (@myloc > c my | oc > BDi , true > y/AB) 


Check patch capacity 


Fin, 


def 


Pb ■ "A2 + — Pb) ■ y/ AtS 


ueciQc wneiner 10 disperse 


/J/12 


def 


£teiVe(g/!(myloc) Pmyloc/ : go£.BC2 


uispersai anempi z 


/JC-2 


def 


con a ^ viv my 10c t.myloc ^ y ..r t , true l> y ./i/j j 


nudiei ui duuii 


Fl 


def 


?/ :FC + (l- ?/ ):0 


Floater survival 


FC Q 


def 


cond (@myloc > c my | oc > FDo,true> y/-Fl) 


Check patch capacity 


FD Q 


def 


pf.FA 1 + (l-p f ):y/.Fl 


Decide whether to disperse 


FA 1 


def 


£teiV«gfc(myloc) Pmyloci' : gol.FC\ 


Dispersal attempt 1 


FC\ 


def 


cond (@myloc > c my | oc >FDi,true[> ^/.F/) 


Check patch capacity 


FDi 


def 


p /: FA 2 + (l-p/) : V-FZ 


Decide whether to disperse 


FA 2 


def 


LfGWeig/i(myloc) /'myloci : goty/ .Fl 


Dispersal attempt 2 



As before, the system can be modeled as the composition of the species as well as the various individuals 
that form the study: 

System^ [(R:[[l]]\ ]J AB:[[£ u l]]\\ ]J Juv.^ulj | ]J Fl:^,!]] . . .)\{re Pu . . .re Pk }. 

l<i<n l b l<i<n [ j l<i'<n} 

Varying the model parameters, e.g. the habitat topology, patch quality and dispersal distance, may allow 
an analysis of the effects of the parameters on patch and metapopulation persistence. 

4 Concluding remarks 

This paper reports on work towards the development of a process-calculus framework for the spatially- 
explicit and individual-based modeling of ecological systems. In related work [1J we have also imple- 
mented a prototype tool and conducted simulations for the spatially-explicit model of [8 ]. In future work 
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we intend to provide optimizations for our tool via an implementation of a spatial extension of the Gille- 
spie simulation algorithm lfl2l I2T1 and by taking advantage of concepts developed in process-algebraic 
frameworks for state-space reduction such as confluence and minimization according to equivalence re- 
lations. At the same time it is our intention to enhance the syntax of PALPS to enable a more succinct 
presentation of systems especially in terms of the multiplicity of individuals. Other possible directions 
for future work include the adoption of continuous time as well as the use of dynamic attributes to allow 
exploring the system while patch quality degrades, temperatures increase, etc. 
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